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We work on the relation between the local thermodynamic instability and 
\ the dynamical instability of large black holes in four-dimensional anti-de Sitter 

^sO ■ space proposed by Gubser and Mitra. We find that all perturbations suppress- 



ing the metric fluctuations at linear order become dynamically unstable when 
black holes lose the local thermodynamic stability. We discuss how dynamical 
instabilities can be explained by the Second Law of Thermodynamics. 



1 Introduction 

Black holes are very interesting objects from their causal structures in general relativity 
to their quantum mechanical properties. To figure out their physical relevance, we need 
to answer if the complete gravitational collapse of a body results in a black hole rather 
than a naked singularity. The conjecture |jj that nature censors naked singularity was 
proposed in this respect. One of motivation of this is from the fact that black holes 
in 4-dimensional asymptotic Minkowski space are stable: linear perturbations around 
black hole solutions do not give any evolution. 

However, it was found in ||, |3| that black strings and p-branes are unstable. The 
basic idea of the Gregory-Laflamme instability is that whatever has the biggest entropy 
is favored. Since a black string has a different topology of horizon as that of a black 
hole and entropy is proportional to the area of horizon, array of black holes has bigger 
entropy than a uncompactified black string when they have the same mass ||]. The 
instability of black strings was shown H, [3| by doing perturbation theory. It is a very 
interesting question to see what would happen during the transition between them. 
It has been argued that violation of cosmic censorship does occur during this process. 
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Recently it has been suggested that a black string settles down to a new static black 
string solution which is not translationally invariant along the string ||. 

Entropy argument used above was revisited by Gubser and Mitra to propose that 
a black brane becomes dynamically unstable when it is locally thermodynamically 
unstable |6|, [/]]. Local thermodynamic stability is defined as having an entropy which is 
concave down as a function of the mass and the conserved charges [p] . This conjecture 
was made from the perspective of AdS/CFT correspondence 0, I0L III], which identifies 



two low energy excitations, both of which are decoupled from supergravity in flat space, 
in two low energy descriptions of superstring theory fI2|| . Some unstable fluctuation 
modes may be excited when there is a thermodynamic instability in the field theory 
and according to AdS / CFT the same thing would happen in AdS J?| . A semi-classical 
proof of above conjecture using the Euclidean path integral approach to quantum 
gravity was given in [|T^] . 

The motivation of Gubser- Mitra (GM) conjecture is that Lorentzian time evolution 
should proceed so as to increase the entropy. In this paper, we find that their argument 
holds for both of the cases in which the metric fluctuations are suppressed and in which 
only the metric fluctuations are turned on. In the former case, any perturbation for all 
equal charges of AdS^-RN solutions becomes dynamically unstable when the system 
loses the thermodynamic stability and all evolutions increase entropy. In the latter 
case, there is no dynamical instability even though the system is thermodynamically 
unstable. The stability can be explained by the fact that entropy would be decreasing 
if the perturbation is unstable. We discuss these in section 2 and section 3. We try to 
explain dynamical instabilities of black holes from the Second Law of Thermodynamics 
in section 4. 



2 Gubser-Mitra analysis and its generalization 
2.1 AdS A -RN black hole 



An electrically charged black hole in the asymptotically AdS^ was found in [14]. Start 



ing from Af = 8 supergravity in 4-dimensions, they gauged the rigid SO (8) symmetry 



of 28 gauge boson |15| and the potential induced from this gauging makes AdS 4 a 
vacuum solution of this theory. The AdS^ black hole solution is made by focusing 
on t/(l) 4 Cartan subgroup of SO(8), which is believed to be a consistent truncation. 
Only three scalar fields of 70 scalar fields in the original theory are kept by working 
in symmetric gauge. The Lagrangian of this truncated theory is, 



2k 2 



where al 



i=l 



A=l 




The metric signature is 



-+) and G 4 = 4. The electrically charged solutions are 
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where the quantities Qa are the physical conserved charges. The mass is JIB 

1 4 



2 + 4 E «a. 
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and the entropy is 



M 



s = nz H JH(z H 



(3) 



(4) 



where z# is the largest root of F(zh) = 0. It is possible to express M directly in terms 
of the entropy and the physical charges in the large black hole limit, M^>L, as 



M 



2tt-2L 2 VS [a=i 



(5) 



We are going to study in the case where all four charges are equal, qA = q- In this case 
the solution can be written in term of a new radial variable, r = z + q, and it becomes 



ds 2 = -fdt 2 + ^+r 2 dn 2 

F * = (6) 



r r 2 L 2 

It is known [|17]] that the consistent S 7 truncation of 11-dimensional supergravity 
is equivalent to M = 8 4-dimensional gauged supergravity Also the equivalence of 
large R-charged black holes in D=4, D=5 and D=7 with spinning near-extreme M2, 



D3 and M5 branes are respectively demonstrated in |TpJ . It is important to check that 
our black holes can be embedded to higher dimensional black objects because GM 
conjecture requests the non-compact translational symmetry. Any instability found in 
the large black hole limit, M^>L in our case implies the instability of M2-brane. 



3 



2.2 Thermodynamic instability and adiabatic evolution 

Local thermodynamic stability is defined as having an entropy which is concave down 
as a function of the extensive variables. This means that the Hessian matrix, 

a 2 s 8 2 S 



H S = 8M 2 8M8Q B ( 7 \ 

M,Qa d 2 S d 2 S K') 

\ 8Qa9M 8Qa9Qb J 

has no positive eigenvalues. It is straightforward to express q a in terms of deriva- 
tives of M(S,Q A ). From this definition, if we introduce the dimensionless variable 
X= — 2 — r for all equal charges, Qa=Q, the thermodynamic instability is present when 

M'SL'S 

X > 1 |0I • We can see that the most positive eigenvetorQ of Hessian increases entropy 
most when entropy is at its extremum. 

S(M + 8M, Q A + 5Q A ) = S{M, Q A ) + \{8M, 5Q a )H^ Qa (8) 

Even though entropy is not at its extremum, we can forget about the first derivative 
parts by energy and charge conservation in microcanonical ensemble. With this it 
was found that in the positive eigenvector direction of Hessian for all equal charges, 
the dynamical instability coincides with the thermodynamic instability with a small 
discrepancy due to numerical errors. It will be interesting to see what would happen 
in other directions. 

Gubser and Mitra analyzed the linear perturbation in which fluctuations of the 
metric are suppressed. The most unstable eigenvector is (SM, 5Qa) = (0, 1, 1, —1, — 1) 
for all equal charges. The condition that the metric decouples at linear order is that 
Qa ' $Qa — 0. In this case 5T a t, vanishes at linear order and we can also see from 
(2) the metric does not change at linear order. It is not difficult to make the linear 
perturbation equations in this decoupling case beyond the eigenvector direction. Our 
original motivation was to see two things: first, even though the system loses the ther- 
modynamic stability, it would not be dynamically unstable if we perturb the system in 
the way of decreasing its entropy. Second, because the eigenvector direction increases 
entropy most, it would be the fastest way of increasing entropy. 

The general perturbation in which the metric decouples for all equal charges is that 

SQ A = (1, a, b, -a - b - l)8Q (9) 

where a, b can be any real numbers. From (2), we can make an ansatz about a relevant 
perturbation 

5(j) 

5<j>i = (1 + a, 1 + b, — a — b) — 

, a, b, — a — b — 1)5F. (10) 

This ansatz relating three scalar fields to one scalar field and four U(l) fields to the 
other U(l) field should be checked if it is consistent with equations of motion and 

lr The positive eigenvector here means an eigenvector with a positive eigenvalue. 
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it turns out to be consistent. The linear perturbation for each is the same up to 
overall factor 



V„V M H 8F 2 

M ^2 A 4 " 



16F^5F^ = (11) 



and the linear perturbation for each Ff$ is also the same up to overall factor 

d5F = d*SF + d8(f)A*F = 0. (12) 

Here F is the background field strength in (6): it is the same for four F^ A \ It is 
remarkable that all directions have the same perturbation equation. The case of a = 1, 
b=—l is the unstable eigenvector and we can see that (11) and (12) are exactly what 
Gubser and Mitra found ||. From their analysis, we can conclude that all perturba- 
tions suppressing metric fluctuation at linear order have dynamical instabilities when 
the system is thermodynamically unstable. We can see that the eigenvector direction 
is the fastest way of increasing entropy not because it evolves fastest but because it 
increases entropy most. It was suspected that some of perturbations (9) would de- 
crease entropy by the continuity of (8) in the case that X is slightly greater than 1, 
but the Second Law of Thermodynamics is not violated in another remarkable waj 
The second derivative parts of (8) with perturbation (9) is 



1 -{8M, SQ a )H s m>Qa L^j = A(M, Q A )( x -l){l + a 2 + b 2 + (l + a + b) 2 }. (13) 

Here A(M, Qa) is a positive when x ~ 1- We can see that all perturbations (9) increase 
entropy when x > !• The factorization of the eigenvalue part X — 1 can explain why 
we have the same linear perturbation equations (11) and (12) for all perturbations (9). 



3 Stability from the metric perturbation analysis 
3.1 Metric perturbation equation 

In the previous section, we observed that there is a dynamical instability in any case 
that the metric perturbation is suppressed. It would be very interesting to see what 
would happen in the case that the metric is also involved. This is very difficult to do 
and we have not succeeded in doing this in the case that all fields are involved. 

In this section, we analyze a simple case: three scalar fields and four U(l) fields 
are suppressed. This perturbation is in the direction of 5M^0, 5Qa=0 for all equal 
charges. From (5) and (8) we can see that entropy is decreasing in this perturbation. 

2 We thank a referee of JHEP for pointing out this and we would like to apologize to the authors 
of H for incorrectly disclaiming their argument in the previous version of this paper. 

The violation of the Second Law of Thermodynamics was suspected from the fact that AdS is not 
globally hyperbolic. We thought the area law for black holes might not hold. However the area law 
needs only a partial Cauchy surface [E0|, so it holds in AdS. cj>i do not satisfy the dominant energy 
condition, but they satisfy the strong energy condition defined in fLs! or the timelike convergence 
condition defined in ^o|. Therefore, the area of horizon can not be decreasing with classical evolutions 
by the area law. 
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Varying (1) yields the equations of motion 

V a V°& + ^ sinh0 4 - 2 £ ar^e^F^) 2 = (14a) 

L A=\ 

d a {^je a ^(F^ ab )) = (14.b) 

R * = - y k{j-i cosh + + 4 £ - ^(f (A) ) 2 } . 

(14c) 

We expect a relevant perturbation is 

5<pi = 

5Fit ] = (15) 

5 g a b = lab- 

We need to check if above ansatz is consistent with equations of motion and it is so if 

7 = ri + l r r = (16) 

where 7 = 7 a a = g ab, jba- The linear perturbation equations from (14. a) and (14. b) are 
automatically satisfied for all equal charges . Now we need to check if (15) is consistent 
with (14. c). The linear perturbation from (14. c) is 



-^V a V c7 - ^Vblac + V b V (c7 a)6 

= Alac + lG( K -Fj'F/ lbd -^ ac F 2 + ^g ac F bd F e dlb ^ (17) 

where we use the totally symmetric notation for ( ), A = — A and F 2 = F ab F ao . Four 
U(l) fields become the same. It is a well-known fact from electro magnetism that if 
there is a source term, a simple gauge choice is not easy. However, if we take the trace 
of (17), it becomes 

- ^V a V a 7 = A 7 - 16 7 F 2 . (18) 

We used the condition, 7^ = — Y r f° r this- Because (18) is a homogeneous equation 
for 7, we can choose a transverse traceless gauge whereby 

V a 7 ab = 

7 = 0. (19) 

See [|Tj| for a detail about this gauge choice. It should be noted that this gauge choice 
is only possible with our ansatz that scalar fiends and U(l) fields are not fluctuating. 
Finally we get the perturbation equation for the metric from (17) following [jlSfl . 

(V b V b + 2A- 8F 2 ) lac - (R c d lad + R a d lcd ) - (2R b J + 32F a b F c d ) lbd = (20) 
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It is the even wave in the canonical form [1191 which is relevant to our case: 

\ 

e iwt P,(cosO). 



( Itt 



lab 



\ 



Itr 
Irt Irr 






(21) 






r 2 k 
r 2 k sin 2 6 / 

It can be easily checked that in this form, r y t t and — Y r have the same equation in (20) 
and this proves that our ansatz (^5|) -(|16D is completely consistent with equations of 
motion. This equality between r f i and — Y r * s expected from the SM perturbation of 
the metric in (6). The equations for j tt and i tr are coupled: 



I r r z sm 6 



7tt + 2fd tltr = 



(22) 



2- + /' 

v r 



d r — 



A2 



f" + 



dg sin Ode 



2f'd t 

Itr H 







/ ' r 2 sin^ J f 2 

Here / is defined in (6) and /' = d r f. Using the form (21) we can make the forth 
order ordinary differential equation for j tr , which is what we are going to study by 



numerics. 



3.2 Numerical analysis 

To carry out a numerical study of (22), we can cast the equation in terms of a dimen- 
sionless radial variable u, a dimensionless charge parameter x, a dimentionless mass 
parameter a, and a dimensionless frequency w introduced in M: 



u 



X 



Q 



, ,2 -r- 1 



a 



w 



wL'i 



(23) 



Then we combine two equations in (22) 



w 2 



f 



f 



T + fd 2 u + \2 J - - f \d u 



al(l + 1 



u 



w 2 



J 



x \7 + fdl+ { 2 u + f ) du 

r ^ X . 2 

J — o h — + w . 




(24) 



+ 2f'w 2 



Itr 



U 

Now we need to specify the boundary condition for ^ rt . We want to place a initial 
data surface touching the horizon at one end and ending on the boundary of AdS at 
the other end. Because AdS does not have a Cauchy surface [p0~| , the future domain of 
dependence of this initial data lies inside of Cauchy horizon of AdS. To define 'small' 
for the perturbation at the horizon, we can use non-singular coordinates, Kruskal 
coordinates [^TJ. Dropping the S 2 piece in (6) and introducing a tortoise coordinate 
r* and Kruskal coordinates (T, R) according to 

dr* 1 
dr J 

e a/'(r*)(±t+r.) = ±T + R. (25) 
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The near-horizon metric is regular 



ds 1 



'-dT 2 + dR' 



(26) 



We can express the Kruskal components 7^ b in terms of the original components 7 a b 

~R 2 ltt - 2fRT ltr + f 2 T\ 
~-TR ltt + 2f(T 2 + R 2 ) ltr - f 2 TB- 
T 2 ltt - 2fRT ltr + f 2 R% 



Itt 



Itr 



Ir 
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f'(r H )(- 


_ T 2 + R 2y 
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f'(r H )(- 


-T 2 + R 2 ) 2 
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f'{r H ){- 


-T 2 + i? 2 ) 2 



(27) 



all of which should be finite as r — > rn on our initial data surface. To avoid the issue 
of mode superposition and a better physical sense in which black holes would form in 
a collapse situation, we would require a surface ending on a future horizon p2| . When 
we approach the future horizon from outside of a black hole region, R = T + 0(r — r#). 
This implies that normalizable wavefunctions jtr must be 0(r - r#) as we approach 
the event horizon. We also want it falling off like 1/r 2 near the boundary of AdS^. 
This boundary condition is taken from the asymptotic behavior 7^ ~ r 3/ y tr . Using 
Maple, we solved (24) numerically. We did not find any unstable mode. At o = 0, 
thermodynamic stability is lost at \ = 1. The smallest w 2 in this case is w 2 = 2. There 
is no normalizable wavefunction with negative w 2 . Negative mode is found at x = 3-7, 
which lies in naked singularity regions and therefore is not relevant. We can conclude 
that in the perturbation, 5M^0, 5Qa=0 there is no dynamical instability, which makes 
sense because this perturbation would decrease entropy if it were unstable. 



4 Discussions and Perspectives 

We have been working on Gubser-Mitra conjecture, which relates thermodynamics to 
dynamics in black holes. When the metric fluctuations are suppressed at linear order, 
all perturbations are dynamically unstable if black holes are thermodynamically unsta- 
ble and all of evolutions increase entropy. However, when only the metric fluctuations 
are turned on, they are dynamically stable. Entropy is not decreasing in this case. 
This result strengthens both the motivation of Gubser-Mitra conjecture, which claims 
that Lorentzian time evolution should go so as to increase entropy, and the validity of 
their conjecture. The Second Law of Thermodynamics is not a rigorous law of nature. 
It is correct in a certain macroscopic time scale. Therefore it is remarkable if dynamics 
obeys this law, even though we are interested in classical instabilities. 

As suggested in [0 , entropic arguments can give a good information not only on the 
existence of dynamical instabilities but also on the direction they point. We can have 
a very heuristic argument why there is a dynamical instability when the perturbation 
(SM, 5Qa) increases entropy. Entropy can be understood as a functional of the fields 
describing black hole. Let Xi(t) denote metric, scalar and U(l) fields. If we write 
down the equation (8) in terms of Aj(t), it becomes 
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S(X i(t0 + *)> = S(X,( to) ) + I (^)^ f f (28) 

Here Xi(t ) is the solution describing our black hole and the derivative of entropy 
with respect to X% is a functional derivative. SXi(t) can be understood as a solution 
of the linear perturbation equation in this case. If we assume that ( sx-sx ) §i ves a 
positive definite inner product between fields, we can see that entropy increases if and 
only if the time dependence of Xi(t) is e lwt with w 2 negative. This argument does 
not hold for a Schwarzschild black hole and this might be due to our assumption on 
the positivity of ( gx-sx- )' It will be a very interesting problem to understand this 
from the quantum theory of gravitation and it might explain the reason why we need 
a non-compact translational symmetry of a black object in Gubser-Mitra conjecture. 

This is also very interesting in mathematical point of view. Suppose we have a 
hypersurface defined by S = S(M,Qa)- Each pertuabation (SM, SQa) corresponds 
to a tangent vector originated from p = (Mq, Qoa) up to a normalization factor in a 
tangent space T p S of the hypersurface at p []. This tangent space can be understood 
as the real projective space 9ftP 4 . The second derivative part of equation (8) is a 
homogeneous polynomial of (SM, SQa) of degree 2 and the zero locus of this polynomial 
gives us an algebraic subvariety of 3?P 4 . This variety separates 9?P 4 into two parts: 
SS > and SS < 0. We can also separate 3?P 4 in another way: A stable region in 
which a perturbation (SM, SQa) gives no dynamical evolution and an unstable region 
in which a perturbation (SM, SQa) gives a dynamical evolution. If our argument in 
the previous paragraph is correct, this means that both of separations are the same 
and it will give a very interesting relation between algebraic equations and differential 
equations^. 

It was argued that a unstable black string settles down to a new static black string 
solution which is not translationally invariant along the string and can be viewed as a 
local entropy maximum but not a global one 0. If the final stage of the evolution is 
a local entropy maximum, we can not say that evolutions from different perturbatoins 
end up with the same final solution. In our case, we have 3 eigenvectors of Hessian with 
the same positive eigenvalues for all equal charge: (a = l,b = —1), (a = —1,6 = 1) 
and (a = —1,6= —1) in (9). Considering the sign of these vectors, there are six most 
unstable perturbation vectors in the tangent space. It would be very interesting to 
see that what will be the final solutions for these perturbations. Finally it is an open 
question to see if evolutions from perturbations around the eigenvectors would result 
in the same final solutions in which the evolutions from the nearby eigenvectors result. 
We leave these questions for future work. 
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